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Abstract. A thin current sheet (TCS), with the width of an order of ther¬ 
mal proton gyroradius, appears a fundamental physical object which plays 
an important role in structuring of major magnetospheric current systems 
(magnetotail, magnetodisk, etc.). The TCSs are nowadays under extensive 
study by means of space missions and theoretical models. We consider a sim¬ 
ple model of the TCS separating two half-spaces occupied by a homogenous 
magnetic held of opposite sign tangential to the TCS; a small normal com¬ 
ponent of the magnetic held is prescribed. An analytical solution for the elec¬ 
tric current and plasma density in the close vicinity of the TCS has been ob¬ 
tained and compared with numerical simulation. The number density and 
the electric current prohles have two maxima each. The characteristic spa¬ 
tial scale zs of the maxima location was investigated as a function of initial 
pitch-angle of an incoming charge particle. The ehect of the thermal disper¬ 
sion of the incoming proton beam have been taken into consideration. 
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1. Introduction 

One of the challenges of modern space physics is understanding of the formation and hne 
structure of current sheets (CSs). The dynamical evolution of CSs and related energetic 
processes are responsible for substorms, solar flares, and other space plasma phenomena. 
In fact, CSs appear an accumulator of magnetic energy, thus the processes of the CSs 
decay consist in the transformation of magnetic energy to kinetic and thermal energy 
of plasma. The experimental information regarding CSs is nowadays provided by space 
missions, specifically those investigating the Earth’s magnetosphere and solar wind (e.g., 
CLUSTER, GEOTAIL, THEMIS etc.). 

The observation of evolution of a CS in Earth’s magnetail shows that during a sub¬ 
storm the thickness of the CS could have a few hundred kilometers [Mitchell, D. G. et. 
ai, 1990; Artemyev, A. V. et.ai, 2010; Nakamura, R. et.al, 2008; Sergeev, V. A. et.al, 
1990], which is of the order of magnitude of a few Larmor radii of proton. Such thin 
current sheets (TCSs) have been also observed in Mercury [Whang, Y. C., 1997], solar 
wind [Gosling, J. T. and Szabo, A., 2008], Jovian magnetodisk [Alexeev, I. I. and Be¬ 
lenkaya, E. S., 2005] and they are expected to play an important role in magnetospheres 
of exoplanets [Khodachenko, M. L. et. ai, 2012]. Numerous in situ observations per¬ 
formed by spacecraft, report the presence of a certain spatial structuring of TCSs, e.g. 
the presence of anisotropic plasma flows and a multi-scale character of the electric current 
distribution nested in the TCSs [Runov, A. et. al., 2006; Runov, A. et.al., 2003]. The 
expected reconnection processes and instabilities occurring in the TCSs are believed to 
be a driver for all the above mentioned large-scale astrophysical energetic and dynamical 
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phenomena. This justihes an actuality of a further investigation of TCSs with a thickness 
of a few thermal proton gyroradii. Analytical methods in that respect appear a powerful 
tool which connects the studied natural phenomenon of a TCS with the backgrounds of 
fundamental plasma physics. As it will be shown in this paper, the estimation of thickness 
of a TCS as a function of the ratio of the incoming flow velocity to the thermal speed can 
be performed with the trajectory method. 

Theoretical studies of a TCS in application to various space plasma problems are related 
with the ideas of Syrovatskii [Syrovatskii, S.I., 1971] regarding the MHD formation and 
dynamics of a TCS near the magnetic held X-line. However, the MHD description is 
inapplicable at small scales comparable to the thermal Larmor radius of a proton. In that 
respect, the plasma kinetic approach is more correct here and may be expected to give 
physically reliable results regarding the processes in the TCSs. 

An attempt of a direct numeric simulation of the motion of particles in a prescribed 
magnetic held geometry of a TCS was made in Eastwood, J. W. [1972]. This approach 
revealed hrst hints of the inhomogeneous distributions of plasma and electric current in 
the vicinity of TCS which may be associated with its hue structure. Later, such numerical 
study approach evolved towards the application of more complex models and algorithms 
such as PIC [Mingalev, O. V. et.ai, 2007] and Vlasov-Maxwell simulation [Schmitz, H. 
and Grauer, R., 2006]. However numerical results not always clearly disclose the physical 
nature of the simulated TCS features. 

The analytical methods of kinetic description of plasma in TCS (developed in parallel 
with the numeric simulations) are related with solving of the nonlinear Vlasov-Maxwell 
equations [Mahajan, S. M., 1989]. In general case, this is a complicated task. At the 
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same time, it is known that the solntion may be expressed via the invariants of motion. 
This fact forms a basis for various analytic approaches to the problem. For example, it 
was shown in Alexeev, I. I. and A. P. Kropotkin [1970] that the magnetic momentum of 
a charged particle is conserved in the case of crossings of an inhnitely thin current sheet. 
This result enabled to calculate the distributions of electric current and plasma in a simple 
prescribed magnetic held topology of a TCS [Alexeev, I. I. and H. V. Malova, 1990]. 

The hrst self-consistent kinetic model of TCS was proposed by Kropotkin, A. P. and 
Domrin, V. I. [1996]. This model used the quasi-adiabatic invariant of proton oscilla¬ 
tion across the TCS {1^ = § V^dz ~ const) without taking into account of the electron 
component. Later on, the development of this approach enabled the inclusion of not only 
the electron component but also strong thermal anisotropy of plasma [Zelenyi, L. M. 
et.al., 2004, 2011; Artemyev, A. and Zelenyi, L., 2013]. It has to be mentioned that 
this self-consistent study of the TCS (based on quasi-adiabatic approximation) requires 
an extensive numerical calculations in order to obtain measurable macroscopic physical 
parameters of the TCS. 

In the present paper we consider an easier analytical approach based on the methods 
applied in Alexeev, I. I. and H. V. Malova [1990]. It gives the consistent results which 
enable a comprehensive characterization of plasma and electric current distributions in 
the TCS without engagement of heavy numerical calculations. The considered simple one¬ 
dimensional case of the magnetic held geometry and the corresponding plasma particles 
motion may be attributed as a zero-order linear approximation of such natural phenomena 
as the current sheets of magnetotail, heliospheric and planetary magnetodiscs. In spite 
of the simplicity, the proposed approach enables characterization of the scaling features 
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of the TCSs which, as it will be shown below, are consistent with the results obtained by 
other researches within more complicated modelling approaches. The presented study is 
based just on the analysis of plasma particles trajectories which are a fundamental issue 
of plasma physics. 

The paper is organized as follows. In Section 2 we describe our analytic solution based 
on the similar approach as in Alexeev, I. I. and H. V. Malova [1990]. In Sections we 
present the comparison of the obtained analytic solution and numeric calculations for the 
density of plasma and electric current. Section4 discusses the issue of the typical scaling 
of the TCS obtained with our analytic solutions and shows its consistency with the known 
results of other authors. The paper hndings are summarized and discussed in Sections. 

2. Analytic solution 

The description of plasma motion in a strong magnetic held, when the scale of the 
magnetic held variation is much larger than that of the quasi-periodic motion of particles, 
is based on the assumption of conservation of the adiabatic invariants. By this, the 
thickness d of a CS may be considered as a scale of variation of the magnetic held, 
whereas the particle Larmor rotation radii (r^’^ for protons and electrons, respectively) 
correspond to scale of quasi-periodic motion. For the case of {r^jf < d} the motion of 
particles can be described by means of the adiabatic theory of a guiding center using the 
decomposition on a small parameter /d < 1. However, in a TCS an opposite condition 
r\ > d usually holds true, and the adiabatic invariants are not conserved [Cary, J.R. 
et.ai, 1986; Sonnerup, B.U.O., 1971], and for the description of particle motions in the 
TCS it is necessary to use a quasi-adiabatic theory [Buchner, J. and Zelenyi, L. M., 1989; 
Zelenyi, L. M. et.ai, 2011]. 
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The approach undertaken in this work is different from the attempts of a self-consistent 
description of a TCS in Kropotkin, A. P. and Domrin, V. I. [1996]; Sitnov, M. I. et.al. 
[2000]; Zelenyi, L. M. et.al. [2004], We study the interaction of moving protons with the 
already existing super-thin current layer (STCL) with e = d/r^ << 1, where d stays for 
the thickness of the STCL and is gyroradius of proton. This super-thin current 

layer is responsible for the existence of the prescribed background magnetic held. The 
goal is to study the distributions of the protons density and electric current in the vicinity 
of the STCL followed from the particle motion analysis. These distributions are related to 
the hne structure of a TCS. The physical nature of the postulated STCL could be related 
with the drifting electrons or with the cold quasi-adiabatic protons. 

The major point of the self-consistent model of TCS developed in Sitnov, M. I. et.al. 
[2000]; Zelenyi, L. M. et.al. [2004] is that the electric current of a TCS and its maximal 
spatial scaling ~ (in case of weak anisotropy) are supposed to be controlled mainly by 
the same sort of particles (e.g. protons). Whereas in our case the presence of a STCL 
with e << 1 created by another sort of particles (e.g., electrons or cold protons) is of 
principal importance. In that sense our approach is to certain extend analogous to the 
’’matreshka” model by Zelenyi, L. M. et.al. [2006]. 

The small parameter e <C 1 related with STCL allows to develop the adiabatic analytic 
approach for the analysis of proton motion. It is used in the asymptotic decomposition 
of physical quantities in the series over the powers of e (i.e. oc e”, where n = 0, 1,2..). 
The vector velocity of proton in a zero order approximation (n = 0) does not change the 
direction when a proton crosses the STCL. Therefore in a zero order approximation the 
STCL can be considered as an inhnitely thin one. 
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The electric field of the charged particles can greatly complicate the analysis of their 
motion even in the zero order of approximation. In that respect we consider a typical 
for the magnetospheric plasma case, when the Debye radius is several times less than the 
Larmor radius of proton (r^ = tl = ^^j^sin6*o, where Oq is the particle pitch-angle; rup 
and e are mass and charge of proton, respectively), so that the effect of the particles own 
electric field is not important and the assumption of a quasi-neutrality of plasma holds 
true. In order to neglect the effects of the particle thermal motion, we consider the case 
when the energy of protons, entering the STCL, is much greater than the thermal energy, 
assuming for proton a simple spiral motion and taking electrons as an incompressible and 
massless liquid. 

In our analysis we use an analog of the standard GSM coordinate system with Z - 
axis perpendicular to the STCL and X - axis directed along the tangential component 
of magnetic field (see in Fig.la). For simplicity we assume the absence of a large-scale 
electric field Ey in the STCL which formally corresponds to the transform into the de 
Hoffman-Teller frame. We also assume that the variation of the normal component of 
magnetic field across the infinitely thin current layer is small AB^/AZ ~ 0. Hence from 
zero divergence of magnetic field the change of tangential component along X has to be 
also small, i.e. AB^/AX ~ 0. Such a configuration is widely used for the analysis of most 
of the astrophysical phenomena, e.g. magnetotail and magnetodisk. 

Within the frame of the made geometry and coordinate system assumptions the compo¬ 
nents of magnetic field are as follows: B^^z) = BxSign{z)] By{z) = 0; Bz{z) = B^ = const, 
so that tan(Q;) = B^/B^ - characterizes the inclination of magnetic field relative STCL. 
To distinguish between the physical quantities in the z > 0 and z < 0 half-spaces we 
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use further on for them the indexes 1 and 2, respectively. We consider a monochromatic 
proton beam which consists of particles having the same pitch-angle 6 and the same en¬ 
ergy S. In this case the whole beam can be described based on the analysis of just one 
proton trajectory. In other words, if the solution for one proton is known, the behavior 
of the beam can be obtained by averaging of this solution over 6 and S. By this, the 
trajectory of a proton in the considered magnetic held geometry has a form of a spiral, 
with the parameters which can be reconstructed from the particle pitch-angle and phase 
of rotation of the proton at the moment before its hrst crossing of the STCL (i.e., the 
boundary between regions 1 and 2). 

After the hrst crossing of the STCL, the proton starts an oscillatory motion across the 
current sheet (i.e., XT-plane), entering consequently in the regions 1 and 2 separated 
by the STCL. During this motion the proton makes a hnite number, n, of crossing of 
the STCL and during this process the particle pitch-angle gradually changes (see Fig.lb). 
Each pair of parameters {0^2^ 4>T2) of oscillating particle for a particular m—th crossing 
(1 < m < n) of the STCL is uniquely dehned by the values of these parameters at the 
previous crossings (up to the very hrst entry into the STCL (m = 1) due to the recurrent 
character of the dependence). The lower index of the phase and pitch-angle (i.e., 1 or 2) 
denotes the region from which the particle approaches the STCL boundary. Note that 
for each crossing of the STCL we distinguish between the phases and pitch-angles for 
the particle entry in the current sheet and escape, 0^2 ^^2) which are marked by an 

appropriate combinations of the indexes (see Fig.lc). 

From the conservation of the velocity vector at the STCL boundary the following rela¬ 
tions between 0™, 0™ and 0™, 0™ hold true: 
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cos(6*^) = — cos(6*|") cos(2a) — sin(6*™) sin(2a) sin((/)™), 


(1) 


sin(0™) cos(0™) = sin(0™) cos(0™), 


( 2 ) 


sm{9^) sin(0™) = cos(6'™) sin(2a) — sin(6'™) cos(2a) sin(0™). (3) 

At the intervals between consequent crossings of the STCL the pitch-angle of the particle 
remains unchanged. By this, the following relation for the particle phase values for the 
consequent crossings of the STCL may be obtained from the analysis of the particle motion 
along its spiral trajectory segment above (or below) the STCL between the crossings 


tan(6*™2 


tan(0™2 


In our analytical treatment we considered a special case of a 1. It corresponds to the 
held topology realized in magnetodisks, as well as in the magnetotail and some regions 
of the magnetopause. The smallness of a enables us to obtain the following approximate 
relations for and from equations (l)-(4) 


= 7r-0j” + 2asin(0j”), 


(5) 


= 27r - 0™ + 2a 


cos(00) 
tan(0™) ’ 


( 6 ) 


= 2ir - - 2 “(j '^0 

tan(6*j”) sm(00') 


(7) 
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Using equations. (5), (6), (7) the change of the particle escape phase and pitch-angle at 
the STCL boundary corresponding to one oscillation period (i.e. between m and m + 2) 
may be calculated as 


( 8 ) 


= 4asin(0f). 


(9) 


The same expressions can be written for the change of the particle entry phase and 
pitch-angle and A0™+2^ 

The ratio of the STCL escape (entry) phase to pitch-angle variations on one oscillation 

^0m+2 

period, i.e. does not depend on a (see equation (8) and (9)). This fact, in view of 

^^ 1,2 

the smallness of a allows to transform this ratio to the differential form: 


dc() 


W) + 


( 10 ) 


dO tan(0) sin(0) 

Note that such transformation is possible only in the case of small jumps of A0™^^ and 
which implies an additional condition on a assumed initially just to be small. This 
condition is a < ^o, where 9q is initial pitch-angle. 

The solution of the equation (10) is 


sin(0) = sin(0o 


^2(71 - 0o) -|- sin(20o) 


( 11 ) 


^^2(71 - 0) + sin(2(/)) 

where (j)Q = and 6q = 61 are the initial escape phase and pitch-angle of the particle 
hrst oscillation cycle across the STCL. The same equation as (11) dehnes also the relation 
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between the particle entry phase and pitch-angle given 00 and 9q stay for their initial 
values. 

The equation (11) can be also found from the simple analysis of the trajectory shape 
and the related areas (full circle or truncated circle) it covers (see Fig.2 a ). In the 
case when areas of the full circle and the truncated circle are equal, the relation be¬ 
tween their radii i?o and Ri can be obtained from simple geometry consideration: Ri = 
Rq\/^/ -y2(7r^^'0)~+”sin^20). On the other hand, taking into account the dependence of 
the particle Larmor radius on the magnetic held and pitch-angle, Ri = m\V\ sm{9)/e\B\, 
we obtain the analogous dependence of sin(6*) on the phase as in equation (11). This 
result may be interpreted as a conservation of the magnetic held hux through particle 
Larmor circle area and then, the area inside a segment of the particle trajectory between 
its consequent crossings of the STCL [Alexeev, I. I. and A. P. Kropotkin, 1970]. 

In the analysis below we will use the following notations: 0 - for the phase of rotation 
of a particle in a homogeneous magnetic held, i.e. on a segment of Larmor spiral between 
crossings of the STCL; 0 - for the escape phase of a particle, crossing the STCL. 

The density of particles with momentum in the range dT = p'^ sin{9)dpd(j)d9 at inhnity, 
is dehned as 

f{p,9,(j))dr = foS{p - po)[6{9 - 9q) + 6{9 - it + 9o)]p‘^ sm{9)dpd(j)d9, (12) 

where /o = no/(47rpQ sin(6'o)), no and po = y/2m£ are density and momentum of parti¬ 
cles, respectively. The phase of rotation 0 and escape phase 0 as the functions of coordi¬ 
nates can be calculated as follows. Expressing 14 on the segment of trajectory {9 = 9{(l))) 
in terms of 0 and 6* as 14 = |l/|(cos(6') sin(a) — sin(6*) sin(0) cos(a)), taking into account 
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the homogeneity of the magnetic held and the dehnition 14 = dz/dt = uidz/dcj), where 
cji is Larmor frequency of proton, one can obtain by integration over 4 from 0 to 4 


z = 


Tl 


sin ( 6*0 


[sin(6*(0)) cos(a)(cos(0) — cos(0)) + cos(6*(0)) sin(a)(0 — 0)]. (13) 


Introducing the projection of Larmor radius on the 2 ;—axis zq = cos(a!), and taking 

into account the smallness of a the expression (13) can be written as 


sin( 0 ( 0 )) ~ 

= ^ 0 —ryryy—[cos(0) — cos(0)J. 


(14) 


sin( 6 >o) 

Substituting the dependence 6{(j)) dehned by equation (11) into equation (14) we obtain 
the relation for 0 to 0 : 


cos( 4 ) = cos( 4 ) — 


z / 2(71 — 0) + sin(20) 


zo y 2(71 — 0o) + sin( 20 o) 
Thus the distribution function (equation (12)) can be rewritten as 


(15) 


/(p,6*,0)dr = foS{p - po)[S{9 - e{^)) + (5(6»-7r + 6»(0))]p^sin(6'(0))^d0dpd6', (16) 

ocp 

where d4>/dcj) follows from (15). 

For integration of equation (16) it is necessary to investigate the boundaries of the 
phase space. In the considered model, two major types of particle trajectory can be 
distinguished (see Fig.l). The hrst trajectory type corresponds to a simple spiral rotation 
in a homogenous magnetic held, before the hrst and after the last crossing of the STCL 
when particle approaches or leaves the STCL. Another type of trajectory corresponds the 
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oscillatory motion of the particle across the STCL accompanied by the particle gradual 
displacement due to the effect of B^. 

Therefore, each particle on its oscillatory part of the trajectory appears to certain 
extend associated with the TCS and contributes the specihc near-by environment which 
may be considered as an extension of the TCS in the range ^ = -^Zq. The range of 
rotational phases of a particle approaching the STCL depends on its z-coordinate. For 
^ > 2zo (i-e. the case of a pure spiral trajectory) the rotational phase changes from 
0 to 271. For z < 2zq the phase of particles changes from 0o(^) to 27r — (1)q{z), where 
0o(-s) = arccos(cos(0o) — z/zq) and 0o is the phase of hrst crossing. By this, 0 changes 
during the particle oscillatory motion across the STCL from 00 tih some maximum value 
(f)m and then back to 0o (Fig.2b). When the 0 reaches its maximum value 0m, the particle 
pitch-angle becomes 7r/2. In view of that one can obtain from equation (11) an equation 
which dehnes the value of 0m: 

2(7r - 0m) + sin(20m) = [2(7r - 0o) + sin(20o)] sin^(6>o). (17) 

Average values of TCS (e.g., density and electric current) in the close vicinity of STCL 
can be calculated as 

<r]{eo,z)>=2 j r]{(p,e)d(p + A J r]{^,e{^),z)d4), (18) 

4>o{z) 0 

where ri{9o, z) corresponds to the particular physical quantity which is averaged. The 
hrst integral in equation (18) corresponds to the simple motion of the proton before hrst 
crossing. The coefhcient 2 before the hrst integral in equation (18) is due to the account 
of both types of protons - the approaching and leaving the STCL. The second integral 
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is related to the protons’ periodical crossings of the STCL, whereas the factor 4 appears 
due to the symmetry of the under-integral function on the intervals (or ■■■ ^/2) 

and 0m---0o (or 7r/2 ... tt — 9q) and account of protons approaching to the TCS from both 
sides (e.g. from the regions 1 and 2). With formula (18) one can calculate the physical 
parameters of TCS for different values of pitch-angle 6q. 

3. Numeric simulation versus analytic solntion 

To check the consistency of the obtained analytic solution we compare it with the results 
of numeric simulations. Taking into account that protons play the major role in dehning 
the parameters of the TCS and assuming their mutually independent motion, we simulate 
numerically the dynamics of protons approaching the STCL from the both sides. 

In case, when the protons move independently from each other the trajectories may 
be calculated independently by numerical method. The standard 4-5 order Runge-Kutta 
numerical method has been used for solving the set of equations of motion. After that 
the calculated vector velocity (for each trajectory) was used to calculate the 9 and 0 for 
the corresponding coordinate . 2 . 

The goal of numerical simulations is to obtain the integral plasma parameters (according 
to equation (18)) dehned with the reconstructed dependence 9 = 9{(j)) (where (j) varies from 
(t^min to (j)max) for £xed Coordinate 2 ;. In course of test simulation runs we found that the 
amount of 30 protons is, in principle, enough to reconstruct 9 = 9{(j)) with a sufficient 
accuracy. To ensure better quality of the simulation and to minimize errors we hnally 
take 60 protons approaching the STCL from the region 1 and another 60 protons coming 
to the STCL from an opposite side, i.e. the region 2. The module of initial velocity was 
assumed to be the same for all particles (i.e. the monochromatic case was considered). 
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To calculate the particle trajectories we used the normalized equation of motion 




[V X B], 


(19) 


In course of the simulation, we use the same coordinate system and the geometry of 
magnetic held as those in the case of analytic treatment and take for the calculations = 
0.01 and Bx{z) = sign{z). The initial velocity of particles was taken V = {V||, Vj_i, Vj_ 2 } 
and |1/| = 1 with the components depending on the initial values 6q, (po as follows 


Vj| = cos(6'o), 


( 20 ) 


V±i = -sin(6'o) cos(0o), 


( 21 ) 


V ±2 = sin(6'o) sin(0o)- (22) 

By this, for the regions 1 (z > 0) and 2 (z < 0) the initial pitch-angle valnes for 
all particles were taken as vr — 6*o and do, respectively. The initial phases for different 
particles approaching the STCL in the same region (i.e. in the region 1 or 2) were taken 
as 00 = 27r/n, where n is nnmber of the particle (n = 1, 2,..., 60). 

The initial position of the particles guiding centers were and —for the particles 
in the regions 1 and 2, respectively. The size of the simulation box in z-direction was 
2) = {-9rL;9ri}. 

Based on all calculated 120 trajectories the distribution function equation (12) was 
reconstructed in the parameter space {6, 0} for the running values of 2 ;—coordinate. Fi- 
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nally the electric current Jy and number density N as the functions of 2 ; are given by the 
integrals: 

(pmax (- 2 ^) 

Jy{z) = Jo J sin^(6*(0)) cos(0)(i0, (23) 

4^min ( 2 ^) 

(pmax (-S) 

N{z) = No j sm{e{(j)))d(j), (24) 

4*min (■ 2 ^) 

where (j)min{z) and (j)max{z) stay for the boundaries of the considered phase-space for 
the corresponding 2 ;. The dimensional coefficients Nq, Jo have the following form No = 
A^o/4vrsin(0o), Jo = eVoiVo where Vq and No are module velocity and the number density 
of incoming flow, respectively. 

To investigate the influence of the particle pitch-angle on the electric current and particle 
number density prohles of TCS along we perform calculations for three different values 
of 60 = {0.2, 0.4, 0.6}. Due to the symmetry of the problem (in the regions 1 and 2), the 
calculation results for only one half-space (e.g., .2 > 0) are presented in Fig. 3. 

Both, the particle number density and the electric current reach the maximum value 
not in the center of the current sheet ( 2 ; = 0), but slightly aside of it (at 2 ; = Zm)- Similar 
effect was discovered also in [Eastwood, J. W., 1972]. The obtained splitting (i.e. shifted 
maximum) of the electric current and particle density prohles is due to the oscillatory 
character of particle motion in the close vicinity {±Zm) of the TCS. Note, that in the case 
of Maxwellian velocity distribution of particles approaching the TCS the effect of splitting 
becomes less pronounced than in the above considered monochromatic case. 

As one can see in Fig.3, the results of the numeric simulation are rather close to the 
analytic solutions given by the equation (18). However, there is some difference in the 
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amplitudes of maxima of the electric current and particle density. That is because in 
obtaining our analytic solution we considered the approximate equations (14) and (15) 
which ignore the terms of the order of a 1. Therefore, the discrepancy between the 
numerical and analytic solutions is of the order of 5 ~ a/tan(6*o) ~ Bz/{ExtoniOo)) 
which yields 5 ~ 0.05 and S ~ 0.025 for 9q = 0.2 and 6*o = 0.4, respectively. Besides of 
that, certain difference between the analytic and numeric solutions is related to different 
account the hrst phase 0o of crossings of protons. The analytic solution was obtained 
for a particular value 0o of the entry phase, the same for all particles, whereas in the 
numeric simulation the particle entry phases are different for different particles. In order 
to compare the numerical and analytical results, an average entry phase for the hrst 
crossing of STCL for the whole ensemble of model particles was used as in fhe analytical 
solution. 

Note, that dimensionless plasma density = N{z > 2rL)/NQ outside the TCS, 
according both the analytic solution and numeric simulation, shows the following de¬ 
pendence on the pitch-angle: 47rsin(0o), whereas the ratio of the TCS density max¬ 
imum Nmax = N{zm)/No to N^o decreases for the increasing Oq = {0.2,0.4,0.6} as 
Nmax/Noo = {1.55,1.25,1.18} (where Zm is coordinate of maximal value of N/Nq). Also 
the ratio of the plasma density Nz=o = A^(0)/Ao in the center of the TSC to N^o decreases: 
Nz=o/Nao = {1.36, 0.93, 0.76} and can even reach the values less than one. Therefore, we 
see that the particles with different pitch-angles contribute differently to the density dis¬ 
tribution in the TCS and its vicinity. Decrease of 6 o results in the compression of the 
TCS, whereas the increase of 6 q produces strong rarehcation of the number density in the 
center (z = 0). 
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Besides of growing Zm, also the value of the electric current maximum increases for the 
increasing pitch-angle: Jmax = = {1.32,1.41,1.45}. Note that the sign of the 

current in the center of the TCS (z = 0) changes for the growing dg: Jz=g = J(0)/Jo = 
(0.6, —0.1, —1.2}. This is because the particles with sufficiently big 6 g create (because of 
their trajectory form) a negative electric current which is opposite to the main current. 

4. The scaling of TCS 

Based on the above presented analytic approach, it is possible to estimate the value of 
the characteristic scale zs = Zml^L of the TCS hue structure (where Zm is coordinate of 
the electric current and plasma density maximum) and its dependence on the pitch-angle 
Qg. We pay attention to the fact that the particles oscillating across the STCL in its close 
vicinity move along the quasi-circular trajectories (see Fig.l), so that the particles with 
higher y— component of velocity give the higher contribution to the electric current. The 
same is true for the particle contribution to the number density. By this, according to the 
equations (2) and (14), Vy is a function of 0, 6*(0) and rotational phase 0. Therefore, the 
location of electric current maximum approximately corresponds 0 = vr and 0(0m) = ii/2. 
In view of that, and assuming for the dehniteness sake 0o = 0, we obtain from the equation 
(14) an expression for the estimation of zs'- 


Zs = 


1 + COs(0m) 


(25) 


sin(6'o) 

where 0m is the solution of equation (17). The Fig.4 illustrates the dependence of zs 
on the pitch-angle 9g as provided by the equation (25). In particular, for the increasing 
9g the coordinate of maximum of the electric current and plasma density in the TCS 
shifts towards larger 2 ; up to the maximum value 2rL. That may be considered as a 
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broadening of the TCS. Note that our analysis is valid under the assumption of a < 6 ^ 0 , 
therefore approaching to 6 *o = 0 corresponds the tangential discontinuity and as result the 
self-thickness became zero. 

In process of crossings the pitch-angle is growing (see. Fig. 2b) and as result the 
maximal Larmor radius {6 = tt/2) can be expressed as sin( 6 *o). In this case the 

expression for the characteristic scale equation (25) yields = 1 + cos(0m), where 

0m is a function of Oq. Taking that the pitch-angle can be expressed as Oq = arctan(Vl/V||) 
and considering the case of a strongly anisotropic plasma beam (i.e., —)■ 0 ) we obtain for 
the TCS scaling the well known relation = (hj_/h||)^^^- This scaling was found 

hrst in Francfort, P. and Pellat, R. [1976] within frame of a similar to our approach. 
Considering the TCS scaling in terms of in the similar case of strong anisotropy one 
can obtain from equation (25) that Zml^L = (Vj_/V||)^/^. The same result was reported 
in Sitnov, M. L et.al. [2000] in the limiting case Bn/B q < Vt/Vr, « 1 (where Id, Vo are 
thermal and flow velocities) which is also supposed to hold true (a < 6 q) in our model for 
case of 00 "C 1 . 

The presented analytic approach enables clear interpretation of the different TCS scaling 
types obtained in Prancfort, P. and Pellat, R. [1976] and Sitnov, M. L et.al. [2000]. In 
particular, the power index 4/3 corresponds the scaling in terms of the TCS internal 
Larmor radius (i.e. r™“^), whereas index 1/3 appears in the scaling dehned in terms of 
the incoming proton Larmor radius {tl). 

Our method enables a generalization for the case when the protons have certain distri¬ 
bution over velocities and pitch-angles {Oq), which is described by the distribution function 
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f{V,9o). The average half thickness of the TCS < 2 ;^ > in this case can be dehned as 
follows 


00 7r/2 

Jdv J deoZmiv,eo)f{v,eo)v^sm{eo) 

< Zm >= - -^- JZ , (26) 

CXD TT/Z 

JdV J dOo f{V,eo)V^sm{eo) 

0 0 

where Zm(y,9) = ^(1 + cos(0m(6'o)))- Note that the dependence of 0m(6*o) is given 

by equation (17), where 00 = 0. For simplicity, we consider here a particular class of 

the distribution function of the following form f{V,9o) = fviV)f0{9o), assuming = 
-{Vp-vf 

e . The last corresponds the shifted Maxwellian distribution with the flow velocity 

Vd and thermal velocity V). Finally, the equation (26) can be rewritten as. 


7r/2 

/ cos(0m(6'o))/0(6'o) sin( 6 'o)d 6 'o 

<z^>=^p,F[-^]{l + Zi - — -}, (27) 

* / /»(^o) sm(fo)d6c 

0 


00 —( 

/C3e V2U 

where pt = — is thermal Larmor radius and the function = -2-, 7 -- 

^V2Vt^ 00 _( yo _(-)2 

/c^e A/2U d( 
0 

can be expressed via the tabulated error integral function {erf{Vu/V^Vt) = 

Vp / \/2Vt 


f exp —C^dQ. The equation (27) describes the dependence of thickness of TCS 
in terms pt as a function of Vd/v^V). 

To investigate the dependence of the scaling of TCS on the incoming proton flow velocity 
Vd we rewrite equation (27) as. 


Zm ^ 
pt 


\/2F[ 


Vp 

\/217 


/C, 


(28) 
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where /C stays for the hnite expression in brackets in the equation (27). The behavior 
of the function F[^^] is presented on hg.5. In this hgure one can see that for large 
values of VD/\/2Vt (> 1.5) the function becomes a linear one (proportional to 

In view of that the scaling expression equation (28) can be used for the estimation 
(with accuracy up to the proportionality coefficient) of energy of the incoming proton 
ffow resulting in particular half thickness of TCS, or vice versa, to calculate the thickness 
of the TCS based on a given energy of the flow, using the dehnitions V) = and 

Vd jnip where 8 is the energy of the flow: 


^ ^ / ‘^8 

(X \l — IC, 


(29) 


Pt V T 

where T is temperature of the incoming proton flow measured in energy units. Specif¬ 
ically, in the case of a proton beam with an energy ~ 12 keV and thermal dispersion 
corresponding T ~ 1.5 keV it follows from equation (29) that oc 4/C. If the dis¬ 

tribution function fg{0) turns /C to ~ 1, then the scale of half thickness of TCS in this 
particular case is ~ 4. Such typical scales are often measured for the TCSs in the terres¬ 
trial magnetotail [Runov, A. et. al, 2006]. 


5. Conclusion 

In this paper we present the results of analytic treatment of proton dynamics in the 
vicinity of TCS and compare them with the corresponding numerical modeling, without 
taking into account of the effect of electrons. Both, numerical and analytic results appear 
in good qualitative and quantitative agreement. That demonstrates the high efficiency 
and potential usefulness of the proposed analytical approach for the analysis of TCSs. 
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It was numerically shown in Harold, J.B. and Chen, J. [1996] that the negative compo¬ 
nent of electric current (diamagnetic current) could play signihcant role in case Vu « Vt- 
The role of dia- and paramagnetic currents was studied theoretically in Zelenyi, L. M. 
et.al. [2004], In particular, in the case of Vd >> Vt, the diamagnetic effects are propor¬ 
tional to Vl and became small due to smallness of Vt. Our results demonstrate similar 
picture, where the negative (diamagnetic) current increases with the increasing pitch- 
angles (i.e. decreasing Vn/Vt) and vice versa. Also the number density in the TCS 
strongly depends on the value of the pitch-angle and the resulting density distribution in 
the TCS is connected with the value of particles pitch-angle, which may (depending on 
particular parameters) look as a density increase or rarehcation. 

Based on the analytic approach it was possible to interpret the obtained distribution of 
the electric current and particle density near the neutral line and to estimate its charac¬ 
teristic scale Zs as a function of particle pitch-angle. The obtained analytic expression for 
zs has been shown to agree in the limiting cases with known estimates of the TCS scaling 
derived under more general assumptions [Francfort, P. and Pellat, R., 1976; Sitnov, M. I. 
et.al, 2000]. Within the frame of the presented analytic approach it was possible to give 
an interpretation of the different TCS scaling types (i.e., the power indexes 4/3 and 1/3) 
which are connected with the scaling in units of differently dehned Larmor radii. Finally, 
the estimations of the average scaling of the half-thickness of TCS, with the inclusion of 
the thermal dispersion of the incoming proton beam gave for the particular (common) 
parameters of the beam the values ~ 4, which are often observed in the Earth’s 

magnetotail. 
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The proposed analytical approach may be efficient in further investigation of TCS hne 
structure in the more general case. In particular, as further development of the pre¬ 
sented study may be construction of self-consistent solutions for the magnetic held and 
electric current by an appropriate adjustment of the proton how and particle pitch-angle 
distribution function. 
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Figure 1. The configuration of the magnetic field (a); A schematic view of the proton trajectory 
oscillating around the TCS with the entry and escape phases denoted for each crossing of the 
TCS (b,c). The lower index of the phase and pitch-angle (i.e., 1 or 2) denotes the region from 
which the particle approaches the TCS, the upper index of phase denotes number of crossing. 


DRAFT 


June 18, 2015, 12:32aiii 


DRAFT 












X - 30 


SASUNOV ET AL.: SCALING OF A THIN CURRENT SHEET 





b) 

Figure 2. The areas of the full and truncated circles covered by the particle trajectory in the 
case of free spiral motion and oscillation across the TCS, respectively a). The behavior of pitch- 
angle 6{(j)) as a function of phase for different values of 6*o b). The red and blue lines correspond 
Oq = 0.2 and 6*o = 0.6, respectively. The phase (pmax is a maximum phase for which particle may 
cross the TCS. 
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0 = 0.2 


9=0.2 



9=0.4 



c) 

9=0.6 




b) 

9=0.4 



d) 

9=0.6 



Figure 3. 


Z/r. 


e) f) 

The distributions of the electric current (a,c,e) and number density (b,d,f) of the 


TCS for the 6q = 0.2, 0.4, 0.6. The red line is result of numerical simulation and the blue line is 
the result of the analytic solution eq. (18). 
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Figure 4. The scaling parameter zs as a function of Oq. 



Figure 5. The dependence of coefficient F on the ratio Vo/v^V^. 
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